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Abstract 

The strong coupling treatment of the Frohlich-type polaronic system, based on a canonical 
transformation and a standard Landau-Pekar type variational wave function, is applied to the 
polaronic system consisting of an impurity in a condensate. Within this approach the Relaxed 
Excited States are retrieved as a typical polaronic feature in the energy spectrum. For these states 
we calculate the corresponding effective mass and the minimal coupling constant required for them 
to occur. The present approach allows to derive approximate expressions for the transition energies 
between different Relaxed Excited States in a much simpler way than with the full Mori-Zwanzig 
approach, and with a good accuracy, which improves with increasing coupling. The transition 
energies obtained here can be used as the spectroscopic fingerprint for the experimental observation 
of Relaxed Excited States of impurities in a condensate. 
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I. INTRODUCTION 



Ultracold gasses have shown to be very useful for experimentally testing predictions from 
many-body theories that arose in the context of condensed matter physics. Their study has 
proven to be especially adequate for strong coupling regimes or when there is strong corre- 
lation pj. Recently it has been demonstrated that if the Bogoliubov approximation is valid 
the system existing of an impurity in a Bose-Einstein condensate (BEC) can be appended 
to this list by a mapping to the Frohlich polaron system 

a, a. 

The Frohlich polaron, 

originally introduced in condensed matter physics, is a quasi-particle consisting of an elec- 
tron with its self-induced polarization cloud. The Frohlich polaron Hamiltonian has resisted 
exact analytical diagonalization and has been submitted to many approximation methods 

for a detailed review of the Frohlich polaron in condensed mater physics we refer to Ref. 

4] and for more specialized topics to Ref. js|). The most successful of these methods is the 
variational path integral approach developed by Feynman in Ref. [6j. A polaronic coupling 
parameter was identified that allows to make a distinction between three regimes: weak cou- 
pling, intermediate coupling and strong coupling. Also the internal excitation structure of 
the polaron was investigated by studying the optical absorption and in the intermediate and 
strong coupling regime a resonance was found sj|- This resonance was associated with 
a transition to a Relaxed Excited State (RES), an excitation in which the induced lattice 
polarization is (self-consistently) adapted to the excited state of the electron. More recently 
the polaron optical absorption was studied numerically with diagrammatic quantum Monte 
Carlo numerical techniques [91. At strong coupling some differences in the linewidth were 

n r 

found between the spectra derived in |8| and [9j. Recently these differences were better un- 
derstood through the introduction of the Extended Memory Function Formalism (see Ref. 



10| ) but they have not yet been completely clarified. Also no material is known that exhibits 



polaronic strong coupling. A better understanding of the intermediate and strong coupling 
regimes is needed, i.a. since polaronic effects may play an important role in unconventional 
pairing mechanisms occurring in high-temperature superconductivity fll |. 

In the present work we start by recalling the polaronic system consisting of an impurity in 
a BEC and some standard results from the strong coupling approximation. Subsequently, the 
standard strong coupling approximation is applied to the system consisting of an impurity 
in a BEC. Within this approximation a study is presented of the Relaxed Excited States and 
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the internal transition frequencies of the BEC-impurity polaron. These results are compared 



with the more detailed study of the internal excitation spectrum of Ref. [12|. Furthermore 
the minimal coupling parameters required for the Relaxed Excited States to appear are 
deduced, and the effective masses corresponding to these states are obtained. 



A. The polaronic system consisting of an impurity in a condensate 

In Ref. [3| it is shown that if the Bogoliubov approximation is valid the Hamiltonian of 
an impurity in a condensate can be written as the sum of a mean field term and the Frohlich 
polaron Hamiltonian H po i, given by: 

% k 

In this expression the first term represents the kinetic energy of the impurity with mass 
mi and position (momentum) operators r (p), the second term is the kinetic energy of the 
Bogoliubov excitations with creation (annihilation) operator at (a£) and the third term 
describes the interaction between the impurity and the Bogoliubov excitations. The Bogoli- 
ubov dispersion is given by: 



% = cWl + (#0A (2) 



where use was made of the healing length of the condensate: £ = 1/ s/WkobbtTq and the speed 
of sound in the condensate: c — h/ (-\/2m B £), with n = N /V the condensate density, cibb 
the boson-boson scattering length and mg the mass of the bosons. The interaction amplitude 
is given by: 



Vk = VNo 



2 



1 1/4 



9 ib j (3) 



where giB is the amplitude of the impurity-boson contact potential which (at low tempera- 
ture) is completely determined by the relative mass m r and the impurity-boson scattering 
length aiB through = 2irh 2 aiB/m r . The Feynman variational path integral technique 
was applied to this specific polaron system in Ref. 13[ and it appeared that in this case the 
polaronic coupling parameter a can be defined as follows: 

a = -% (4) 
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As in the case of the solid-state Frohlich polaron a transition was found between a weak 
coupling regime and a strong coupling regime. Furthermore it was noted that the polaronic 
coupling parameter (Hj) can be tuned externally by a magnetic field through a Feshbach 
resonance (see e.g. |14j). With this technique a coupling parameter a of the order of 10 is 
expected to be experimentally feasible. This suggests that by examining an impurity in a 
condensate it may be possible to experimentally reveal the internal excitation structure of 
the polaronic strong coupling regime for the first time. 



B. Standard polaron strong coupling treatment 

In this section a summary of the polaronic strong coupling approximation is presented 
which was introduced by Landau and Pekar for the description of the ground state in Ref. 



Id. These r esu.ts were shown to he accurate in the asymptotic strong coupling limit for 
a symmetrical, exactly soluble one dimensional polaron model in Ref. [16[. The strong 
coupling formalism was later extended for the calculation of the energy level of the lowest 



Relaxed Excited State in Refs. 



17 



18|. 



The product Ansatz is assumed which means that the wave function of the system \$) 
can be written as the product of the wave functions of the impurity and the Bogoliubov 
excitations \<p): = \<p). The wave function |0) is assumed to be of the form: \<p) = 
U |0), with |0) the vacuum and U a canonical transformation: 



U = exp 



(5) 



with {fg-} variational parameters. If the expectation value of the polaron Hamiltonian with 
respect to \^f) is minimized as a function of the variational functions {f^} the following 
standard expression is found for the energy: 



Vr 



E = K-Y,— 



\Pk\ 



(6) 



where the kinetic energy K and the Fourier transform of the density of the impurity were 
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introduced: 



Pk = 



2mi 

ik.f 



-07 
$1 



(7) 
(8) 



Expression (jfJJ) is strictly speaking an upper bound for the groud state energy. When con- 
sidering excitations of the impurity's wave function, it is assumed that expression <^j, with 
appropriate pr, gives a good approximation for the energy of the excited polaron. 



C. Effective mass 

The above formalism can be extended to allow for a calculation of the strong-coupling 
polaron effective mass 
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20] . The total momentum V of the polaron system, given by 



V =p + Ylk ^otag, commutes with the polaron Hamiltonian (JT]) and is conserved. After 
introducing a Lagrange multiplier v the following expression has to be minimized: 



H (v) = H pol - v. 



(9) 



V 



with v a Lagrange multiplier which physically corresponds to the averaged velocity of the 
polaron. The polaron effective mass m* is then determined from the relation between the 
momentum and the averaged velocity of the polaron: V = m*v. The trial wave function for 
a system in which the polaron moves with velocity v is, i.e.: 



exp 



tmjv.r 



\1>i) 10) 



(10) 



If the expectation value of H (v) is taken with respect to |^ > and then minimized as a 
function of the variational functions {f^\ the following expression is found for the energy: 



E(v) 



K 



nijv 2 , ^ \- |Vfc| 2 |p^| 2 



+ V.V-J2 



hbJt — hv. k 

k k 



(11) 



The next step is the minimization with respect to the Lagrange multiplier v. If a Taylor 
expansion is performed for small v an expression is found for the effective mass: 



mi 



\ v k\ \Pk 



(12) 



with i,j = x,y,z. Expression ( TT2|) results in an anisotropic effective mass for a general 
density pg. Since there are no preferred directions in the system, the physically realized 
effective mass is the lowest eigenvalue of the tensor m*-, since it results in the lowest energy. 



D. Variational wave function for the impurity 

In the present strong coupling approximation it is assumed that the impurity is localized. 
The self-induced potential is approximated as quadratic which means that the impurity wave 
function is chosen to be an eigenstate of the harmonic oscillator: 

VW.(A,*0= II (2^!)- 1/2 ^l/^ (~) exp (~), (13) 

with H ni the Hermite polynomials, Hi — 0, 1, 2, ... and A the oscillator length which is taken 
as a variational parameter. Physically the oscillator length A can be understood as a measure 
of the extension of the wave function, a small A corresponding to a small wave function and 
thus a highly localized state. The kinetic energy flTj) is given by: 

*»< a >4(H^' (14) 

with n = n x + n y + n z . The Fourier transform of the density ([H]) can also be calculated as a 
function of A and k and it turns out to be a function of q = Xk, i.e. (A) = p(q). 



II. STRONG COUPLING APPROXIMATION APPLIED TO AN IMPURITY IN 
A CONDENSATE 

In this section we introduce the Bogoliubov dispersion ()2]) and the interaction amplitude 
03]) for the polaronic system consisting of an impurity in a condensate. If the impurity wave 
function is assumed to be a harmonic oscillator eigenstate ( TT3|) the following expression is 
found for the energy ([6]) (with polaronic units, i.e. h = £ = mj = 1): 

E^) = ^A- + A-—Jdk^ (15) 
v ; 2A 2 V2 j 2i 2 j 2 + P' K J 

where a is the polaronic coupling parameter (jl]) and fi is given as: 

{l + m B f 

p= — . 16 

4m R 



From (|T5|) . it can be seen that the effective coupling parameter is actually a product of a 
and the mass imbalance ft between bosonic atoms and the impurity atoms. The effective 
coupling constant afi can be enhanced by choosing strongly unequal masses, making fi large. 
In what follows, we study the system properties as a function of a/i rather than a. Since 
the density of the harmonic oscillator eigenfunctions is a function of Xk only, (j!5p can be 
rewritten as: 

™-MH-&J*£¥? (17) 

Minimization of ( 1T7|) results in the following condition on the variational parameter A: 

This equation determines a minimal value for a fx below which no solution is found, i.e. the 
requirement to find a solution is: 



au > min 

\ell+ 



7T 2 (n + 3/2) 



(19) 



. 2A3 (/^S? + A/*3SS*)j 

Values of afi not satisfying (fl9l would give A — > oo and thus a completely delocalized 
impurity meaning that the self-induced trapping potential is too shallow to have this type of 
wave function as a bound state. For the ground state this behavior was already observed in 
Ref. [lj where above a critical a a state with a large effective mass and a localized impurity 
wave function was found. 

Introducing the Bogoliubov dispersion (j2J) and the interaction amplitude (j3J) in the effec- 
tive mass (TT2]) results in (in polaronic units): 



Aa fim 2 B \ f p 2 (q) q iqj 
vr 2 J ^(2A 2 + ? 2 ) 2 ? 2 



K = i + / (20) 



III. RELAXED EXCITED STATES 

If an excited state wave function is considered for the impurity, the resulting polaronic 
state is a Relaxed Excited State (RES). This corresponds to an excitation of the impurity 
in the self-induced potential adapted to the excited state's wave function. In this section 
the numerical solutions for the strong-coupling approach are examined for the ground state 
and the three lowest Relaxed Excited States. Polaronic units are used, i.e. h = mi = £ = 1. 



7 



10 20 30 40 50 

FIG. 1: The variationally determined energy levels in the self-induced potential as a function of 
the coupling parameter a[i for the ground state ("0"), for the first Relaxed Excited State ("1") 
and for the two types of second Relaxed Excited States ("2a" and "26"). 
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FIG. 2: The variationally determined harmonic oscillator length of the impurity wave function as 
a function of afx for the ground state ("0") and the three lowest Relaxed Excited States ("1", "2a" 
and "26"). 

For the ground state, n = has to be considered. The first RES corresponds to n — 1, 
i.e. one of the n^'s is taken to be 1. The second RES, corresponding to n = 2, can be realized 
in two distinct ways: either one rij equals 2 or two n^s equal 1 which will be indicated as 
2a and 26, respectively. 

In figure [1] the variationally determined energy levels are shown as a function of a.\x for 
all states discussed above and in figure |2] the variational optimal oscillator length of the 




8 



ground state 
first RES 
second RES (2a) 
second RES (2b) 



(an) 



exist 



A 



(afi) 



stable 



A 



3.83302 0.8213 

9.10277 0.7646 

15.3764 0.7400 

18.0197 0.6852 



3.56641 1.2668 
3.55492 1.1850 
14.4869 1.1495 
17.1152 1.0570 

TABLE I: The values of afi above which the self-induced potential supports the indicated levels 
is given, along with the corresponding values of the variational parameter A (cf. eq. (19)). The 
last two columns give the values of afi (and the corresponding A) above which the indicated levels 
become stable. 

impurity is shown. Note that the energies of the excited states possess no variational upper 
bound character. To the left of the curves, the coupling afi is too weak for the self-induced 
potential to support these levels. For couplings just large enough for the level to exist, it is 
metastable. There are, for each state, two special values of afi'- one at which the self-induced 
potential allows for that state to exist, and one at which the polaron energy (jBJ) becomes 
negative, leading to a stable state. We list these values together with the corresponding 
optimal value of the variational parameter A in table [H From figure [2] it is clear that an 
increase of afi results in a decrease of A, corresponding to a stronger self-induced potential 
and thus to a more localized impurity within the polaron. Furthermore we see that at a 
given afi different oscillator lengths are found for different states which means that not only 
the wave function of the BEC relaxes but also the impurity's wave function is adapted. 

In figure [31 (m* — l)/m 2 B is shown, where m* is the effective mass for the states under 
consideration. The largest effective mass is found for the ground state and it is observed 
that the effective masses corresponding to the two possibilities for the second RES cross 
at a coupling parameter afi ~ 28.7. In figure H] the transition energies are shown for the 
transition from the ground state to the Relaxed Excites States under consideration and from 
the first RES to the two types of the second RES. For larger afi more Relaxed Excited States 
are expected to exist and as a result there will be more possible transition energies. 

Recently a more detailed study of the excitation structure of the polaron system consisting 
of an impurity in a condensate was performed by calculating the response of the system to 
Bragg spectroscopy within the Mori-Zwanzig projection operator formalism 12j. In this case 
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FIG. 3: This figure shows (m* — l)/m 2 B where the effective mass m* is given by eq. (20) for the 
ground state ("0") and for the three lowest Relaxed Excited States ("1", "2a" and "26"). The 
values of (ctfi) at which the self-induced potential is deep enough to support these states are also 
indicated (cf. Tabled]) by dotted vertical lines. 
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FIG. 4: Transition energies from the ground state ("0") to the Relaxed Excited States ("1", "2a" 
and "26") and from the first RES ("1") to the two types of the second RES ("2a" and "26"). The 
critical coupling parameters {o-n) exist from tableUfor the Relaxed Excited States are also indicated 
by dotted vertical lines. 



the system consisting of a lithium impurity in a sodium condensate was examined which 
results for the present treatment in fi = 1.5209 for the mass parameter f lT6|) . With this value 
of n the present theory is able to describe the strong coupling ground state from a ^ 2.5. 
The self-inducedpotential becomes deep enough to support a first Relaxed Excited State at 



«Ri6. In Ref. 



121 ] the transition to the first Relaxed Excited State is already seen at a = 4 
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and for larger a a larger transition frequency is found, a behavior which is also observed 
in figure HI This shows that the lowest a at which the Relaxed Excited State appears is 
slightly smaller in the treatment of Ref. [12j and somewhat larger in the present strong 
coupling treatment. A comparison of the transition frequency from the ground state to the 
first Relaxed Excited State can be made for a larger than 6. At a = 8 this results in a 
difference of about 10% and at a = 10 a deviation of about 5% is found. This trend of a 
better correspondence between the present strong coupling approximation and the arbitrary 
coupling treatment of Ref. [12] at higher coupling parameter is consistent since the present 
strong coupling theory is expected to be more accurate at larger a. 

In the present strong coupling calculation the effect of the mass imbalance is contained 
solely in the factor /j which grows as the masses differ more, leading to a larger a/x. This 
observation provides an additional tool for facilitating the probing of the strong coupling 
regime, namely choosing the masses as different as possible: a lithium impurity in a rubidium 
condensate will correspond to a more strongly interacting system at the same (Feshbach 
adapted) scattering lengths than a lithium impurity in a sodium condensate. 

IV. CONCLUSIONS 

In this work the standard strong coupling approach to polaron theory is applied to in- 
vestigate the strong coupling regime of the polaronic system consisting of an impurity in a 
Bose-Einstein condensate. Within this formalism the critical coupling parameters required 
for Relaxed Excited States to appear were deduced, and the effective masses of these states 
were calculated. For impurity atoms in a condensate, these states can be experimentally 
most easily probed by laser spectroscopy. The calculation of the transition energies between 
different Relaxed Excited States presented here offers a straightforward way to estimate 
the excitation spectrum. Comparison with the path-integral results for the spectral density 



obtained within the Mori-Zwanzig framework (see Ref. 12j) shows that the current method 
quickly becomes very accurate (<5% for a > 10) as the coupling strength grows. Further- 
more it was shown that the polaronic strong coupling regime is easier to reach when the 
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masses of the bosons and the impurity differ more from each other. 
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